Draft version February 1, 2008 

Preprint typeset using LAT^X style emulatcapj v. 04/03/99 



HYDRODYNAMICS OF THE STREAM-DISK IMPACT IN INTERACTING BINARIES 

Philip J. Armitage 1 and Mario Livio 

Space Telescope Science Institute, 3700 San Martin Drive, Baltimore, MD 21218; 
armitage@cita.utoronto.ca, mlivio@stsci.edu 
Draft version February 1, 2008 



Or 
C/3 



> 

in 
m 
O 

a\ ■ 

o : 
r» . 
o\ ■ 

^: 

Or 

6 ■ 



X 



ABSTRACT 

We use hydrodynamic simulations to provide quantitative estimates of the effects of the impact of 
the accretion stream on disks in interacting binaries. For low accretion rates, efficient radiative cooling 
of the hotspot region can occur, and the primary consequence of the stream impact is stream overflow 
toward smaller disk radii. The stream is well described by a ballistic trajectory, but larger masses of 
gas are swept up and overflow at smaller, but still highly supersonic, velocities. If cooling is inefficient, 
overflow still occurs, but there is no coherent stream inward of the disk rim. Qualitatively, the resulting 
structure appears as a bulge extending downstream along the disk rim. We calculate the mass fraction 
and velocity of the overflowing component as a function of the important system parameters, and discuss 
the implications of the results for X-ray observations and doppler tomography of cataclysmic variables, 
low-mass X-ray binaries and supersoft X-ray sources. 

Subject headings: accretion, accretion disks — instabilities — hydrodynamics — novae, cataclysmic 
variables — binaries: general 



1. INTRODUCTION 

Observations of cataclysmic variables, low-mass X-ray 
binaries, and supersoft X-ray sources show that accretion 
disks in these systems display pronounced deviations from 
axisymmetry, for which the most obvious agent is the im- 
pact of the gas stream from the companion star onto the 
disk (for a review, see e.g. Livio 1993). The complexity of 
the stream-disk interaction mandates the use of hydrody- 
namic simulations to fully explore its consequences. 

The dynamics of the accretion stream up to the point 
of impact with the disk was described by Lubow & Shu 
(1975, 1976). The gas stream leaves the inner Lagrange 
point in approximate vertical hydrostatic equilibrium, and 
flows inward along a ballistic trajectory. As it does so, in- 
ertia precludes the stream from adjusting fast enough to 
maintain hydrostatic balance, so that upon reaching the 
disk edge the stream is substantially more extended ver- 
tically than the local hydrostatic value. Gas at several 
scale heights above the disk midplane may then be able 
to overflow the disk rim towards smaller disk radii. Some 
possible observational signatures of this were discussed for 
X-ray binaries by Frank, King & Lasota (1987), and for 
cataclysmic variables by Lubow (1989). 

In this investigation we extend the study of the stream- 
disk interaction by performing three dimensional calcula- 
tions of the stream-disk interaction using the ZEUS hy- 
drodynamics code (Stone & Norman 1992a, b). Our goals 
are to quantify the degree of stream overflow, to explore 
the structure of the stream-disk impact region at the disk 
rim, and to investigate the efficiency of entrainment once 
the gas stream is overflowing the disk. These theoreti- 
cal questions correspond closely to observational issues; 
how much absorption can be generated by an overflowing 
stream, how extended is the hotspot region, and can the 
stream-disk impact generate significant emission at veloc- 
ities intermediate between those of the disk and a ballistic 
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gas stream? 

The layout of this paper is as follows. In §2 we dis- 
cuss and present the results of hydrodynamic simulations 
of the stream-disk impact. We compare our results with 
our previous calculations which used smooth particle hy- 
drodynamics (Armitage & Livio 1996), and with two- 
dimensional calculations (Rozcyczka & Schwarzenberg- 
Czerny 1987). §3 quantifies the degree of stream over- 
flow as a function of the parameters of the stream and 
disk, while §4 analyzes the calculated velocity profile of 
the stream and entrained disk gas. §5 discusses how the 
simulations compare to current observational results, while 
§6 summarizes our main conclusions. 

2. HYDRODYNAMIC SIMULATIONS OF THE INTERACTION 

Hydrodynamic calculations of the interaction of the 
stream with the accretion disk were carried out with the 
ZEUS-3D code developed by the Laboratory for Compu- 
tational Astrophysics (Clarke, Norman & Fiedler 1994). 
ZEUS-3D is a finite difference code utilizing an artificial 
viscosity to capture shocks, which has been extensively 
used and tested for astrophysical fluid dynamics problems. 
Details of the algorithms used in the code are described 
by Stone & Norman (1992a, 1992b). 

2.1. Initial and boundary conditions 

The simulations are performed in Cartesian [x, y, z) ge- 
ometry, with cubic cells. The computational box sur- 
rounds the impact point of the stream and the disk, and 
is aligned with the y axis in the radial direction and the 
z axis aligned vertically. Symmetry about the disk mid- 
plane is assumed. The inflow of the stream and disk are 
prescribed as boundary conditions, and the calculations 
are run until the mass and vertical kinetic energy within 
the box reach an approximate steady state. This requires 
~ 2 flow crossing times of the grid. Gravity from the ac- 
creting star only is included. The grid covers 50 % of the 
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radius of the disk, and about 45° in azimuth. 

The disk prior to impact with the stream is taken to be 
in hydrostatic equilibrium in the vertical direction, which 
together with an isothermal vertical structure corresponds 
to a gaussian density profile with scale height H^. At the 
outer edge of the disk, i? ut, the sound speed, c s is set 
such that the Mach number of the flow is 30. The disk 
is always taken to be isothermal in the vertical direction, 
for most calculations we additionally assume that the disk 
is isothermal radially, so that -ffd(-R) oc i? 3 / 2 . The flar- 
ing of the disk implied by this means that the strongest 
interaction with the overflowing gas will be near the disk 
rim. 

Since the hydrodynamics of the impact are essentially 
those of hypersonic flow past an obstacle, the sharpness 
of the radial truncation of the disk at the outer edge may 
be expected to be of some importance. Tidal torques, and 
possibly the ram pressure of the stream itself, will tend to 
produce a sharp edge to the disk, but this will be opposed 
by the turbulent motions giving rise to the viscosity. Con- 
sequently, we assume a gaussian fall-off of the disk mid- 
plane density pdo beyond R ou t, with a scale length equal to 
the vertical scale height H^. Interior to i? ut the midplane 
density is constant with radius. 

The initial conditions for the gas stream are similar to 
the results of the integration of the stream dynamics de- 
scribed by Lubow (1989). We assume that the stream 
density falls off with distance from the stream center as 
a gaussian, with scale length H s , and central density p s o. 
The initial vertical velocity is taken as v z oc —z, with the 
constant of proportionality equal to 0.6c s , implying mod- 
est convergence of streamlines towards the disk midplane. 
The Mach number of the stream at the disk edge is 30, 
and the angle between the stream and the radial direction 
at the impact point is 15°. The stream gas is taken to 
have the same temperature as the outer edge of the disk, 
for many systems this will be a reasonable assumption. In 
cases where the outer parts of the disk are much hotter 
than the surface layers of the secondary, ballistic stream 
overflow would of course be inhibited. 

Slow inflow through a viscous disk implies that pdo 3* 
p S Q. To within factors of order unity, a steady-state disk 
described by the Shakura-Sunyaev (1973) a-prescription 
for the viscosity has, 



Pdo \v s J \H S 



(1) 



with v s the stream velocity at the impact point. For the 
parameters of our calculations, and assuming a typical a 
value of 0.1 for interacting binary disks, p s o/pdo ~ 10~ 2 , 
and we use this density ratio in the simulations. 

Reflecting boundary conditions are imposed at the disk 
mid-plane, z = 0. Inflow boundary conditions are speci- 
fied where the disk flows into the computational volume, at 
x = 0, and at y = over the area where the stream enters. 
Over the remainder of the y = plane, and for the other 
three faces of the computational box, outflow boundary 
conditions are used. Additionally, for the adiabatic sim- 
ulations, we allow outflow through the x = face where 
the calculated inflow density for the disk is below 10~ 12 
of the central disk value (i.e. at high z). This is required 
in order to prevent artefacts arising near the boundary at 



several scale heights above the disk midplane. 

For reference, the calculations use second order (van 
Leer 1977) interpolation for all advected quantities, a di- 
mensionless coefficient of artificial viscosity C2 = 2.0, and 
a Courant number of 0.5. 

2.2. Cooling 

We have run simulations with three equations of state; 
isothermal, adiabatic (with no cooling), and one with opti- 
cally thin radiative cooling included. The isothermal and 
adiabatic with no cooling cases bracket the range of pos- 
sible behavior, the calculation including cooling is more 
realistic although there are many complexities that we do 
not attempt to model in these calculations. 

If cooling is so rapid that the equation of state is effec- 
tively isothermal, then the main parameters of the problem 
as set up here are the Mach number of the flow, the ratio 
of disk to stream density, and the ratio of disk to stream 
scale heights. The results in this limit are independent 
of the absolute disk density, or equivalently of the mass 
transfer rate in the stream. This degeneracy is broken 
when more realistic cooling is included, dense gas then 
cools faster than diffuse material, and the results depend 
on the actual value of M. 

For a simple description (following Blondin, Richards & 
Malinowski 1995) the cooling time is given approximately 

by, 



tc. 



™A(T gas ) ' 
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where ks is the Boltzmann constant, n is the number den- 
sity of electrons, and A is the cooling rate in ergs cm~ 3 
s . We use for A(T gas ) the cooling function given by Dal- 
garno & McCray (1972, see also Cox & Daltabuit 1971), 
for which t c (defined here as the time required to cool to 
7gas/2) is given at a fiducial temperature of 10 6 K by, 
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(3) 



The qualitative hydrodynamics of the flow will then be de- 
termined (Blondin, Richards & Malinowski 1995) by the 
ratio of t c to the dynamical time at the edge of the disk, 
td ~ fi -1 , where 51 is the Keplerian angular velocity. For 
many systems, td will be of order 10 3 s or longer. 

For an estimate, we assume that the emission from the 
hotspot arises from material at comparable density to that 
in the stream (the post shock density will be higher than 
this, so cooling will be more rapid in the optically thin 
limit). Then to order of magnitude, 



n ~ 2 x 10 
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Comparison of this equation with equation (3) indicates 
that in the regime where optically thin cooling described 
by a cooling function is a valid approximation, typical ac- 
cretion rates imply a cooling time t c -C i<j. Thus, at low 
accretion rates we expect the gross hydrodynamics of the 
stream-disk impact to be well described by an isothermal 
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equation of state. The flow may be expected to depart 
from isothermal behavior in the low density regions sev- 
eral scale heights from the midplane, and to explore this 
regime we have run several simulations including optically 
thin cooling. 

2.3. Dependence of the effective equation of state on the 
accretion rate 

The above argument suggests that the accretion rate 
will generally be high enough that cooling will be rapid, 
provided that the emission is optically thin. In this regime 
an isothermal equation of state should provide a good de- 
scription of at least the higher density regions of the flow, 
within a few scale heights of the disk midplane. At higher 
accretion rates the assumption that the optical depth to 
the emitting regions of the hotspot is small will break 
down, and qualitatively we expect that the hydrodynam- 
ics will be better approximated by an adiabatic equation 
of state. As we demonstrate later, the flow in these two 
regimes is very different - here we estimate where the tran- 
sition between them occurs. 

The radiation emitted as the shocked stream gas cools 
must escape, at a minimum, through the hot (~ 10 6 K) 
shocked layer and the relatively cool (~ 10 4 K) inflowing 
stream material. We consider first the opacity due to the 
cooler stream. The temperature of the stream at the outer 
edge of the disk will be comparable to the surface temper- 
ature of the mass losing star, with a modest enhancement 
from compressional heating as the stream converges on the 
disk midplane (in some systems irradiation of the stream 
by the disk or hotspot might also be important - this would 
be particularly significant if the stream gas was able to be 
highly ionized prior to striking the disk). The dominant 
opacity source is therefore likely to be H~ scattering. An 
analytic approximation to the opacity in this temperature 
range is given by Bell & Lin (1994), in units of cm 2 g _1 , 

k = «V /3 T 10 (5) 

with Ko = 10~ 36 . The column density for escape of the 
hotspot radiation will evidently be ~ p s oH s , for which the 
optical depth is, 



Imposing the condition that t = 1 then gives an estimate 
for the critical accretion rate for the transition between 
isothermal and adiabatic behavior. Using the numerical 
values adopted previously we obtain, 

M crit ~ 2 x 1(T 9 ( - H * ) ( ^—r] (7) 

V5 x 10 9 cmy V300kms-V 

For accretion rates below this value cooling should be ef- 
ficient, above it the radiation will be trapped by the in- 
flowing stream gas and the initial cooling of the hot spot 
region will occur via adiabatic expansion. 

The opacity of the hot shocked layer can also limit the 
applicability of the efficient cooling assumption. At M cr it, 
the central stream density for the parameters used pre- 
viously is ~ 5 x 10~ n gcm~ 3 . While the shock is still 



isothermal, the compression could be as large as the Mach 
number (~ 30) squared, giving a shocked layer with a den- 
sity of <~ 5 x 1CP 8 gcm~ 3 . For an opacity due to electron 
scattering, k = 0.348, and r = 1 requires a thickness of 
the hot layer of order 10 8 cm. As this is indeed roughly 
the thickness expected of the shocked region, it suggests 
that the critical accretion rate for the emission to become 
optically thick from the opacity of the shocked layer itself 
is comparable to, or somewhat less than, the critical value 
estimated above. We conclude that, with considerable un- 
certainties, both arguments give M C rit ~ 10~ 9 M yr _1 . 

Clearly the above analysis is extremely crude, and no 
substitute for a proper radiative transfer calculation of 
the hotspot region. A priori it would definitely be wise 
to keep an open mind as to whether a system with a given 
accretion rate would be better described theoretically with 
isothermal or adiabatic models. However it suggests that 
the hotspot region in low accretion rate cataclysmic vari- 
ables might well be able to cool efficiently, whereas nova- 
like variables and supersoft X-ray sources with accretion 
rates much larger than M C rit are almost certainly unable 
to do so. 

In the following section we present models for both lim- 
iting cases. For low accretion rate systems we use an equa- 
tion of state that is either isothermal, or which includes 
radiative cooling in the optically thin limit. These are 
qualitatively very similar, and differ mainly in the very 
low density regions where optically thin cooling becomes 
ineffectual. For models intended to represent better higher 
accretion rate systems we use an adiabatic equation of 
state without any cooling. Real systems, where cooling 
will occur in the low density parts of the flow, would be 
expected to have intermediate properties between the lim- 
iting extremes presented here. 



2.4. Structure of the impact region for isothermal and 
adiabatic equations of state 

Figure 1 (Plate XXX) illustrates the primary differences 
between the simulations using isothermal (or efficient ra- 
diative cooling) and adiabatic equations of state. The left 
panels show isodensity surfaces for an isothermal calcula- 
tion, plotted at a density that is 10~ 3 of the central disk 
density, or 1/10 of the central stream density. The colors 
represent the radial velocity on this surface. 

For this equation of state, the uppermost part of the 
stream where the density significantly exceeds that of the 
disk is able to flow smoothly over the rim of the disk in 
the manner envisaged by Lubow & Shu (1976), and Frank, 
King & Lasota (1987). We do not observe any instabili- 
ties that lead to rapid entrainment of the stream as it 
flows over the disk surface, though the leading edge of the 
stream is deflected by the ram pressure of the disk gas. 
More interestingly, there is a fairly broad region of the 
disk, downstream of the region where the stream is over- 
flowing, which possesses significant inward velocity. This 
is material that has been brought to an intermediate ve- 
locity between that of the stream and the disk, as a result 
of the stream-disk interaction. In our simulations includ- 
ing radiative cooling this gas is hot, and thus a potential 
source of emission with a radial velocity that is highly su- 
personic, but less than that of the ballistic stream. 
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The right panels of the Figure show the results for an 
adiabatic equation of state. Here, the isodensity surface is 
plotted at a density that is 1 CP 2 75 of the central disk den- 
sity, and the colors represent log(c s ) on that surface. The 
details of this calculation will be discussed in Section 2.6, 
but it is obvious from the Figure that in this case the upper 
regions of the stream are not able to flow unimpeded over 
the surface of the disk. Rather, the qualitative impression 
is that the impact leads to a hot bulge extending along the 
rim of the disk. This material is not in hydrostatic equi- 
librium, and expands outwards downstream of the impact 
point. Material is still overflowing towards smaller radii, 
but it is not in the form of a coherent stream. We also note 
that even in the absence of radiative cooling, the temper- 
ature on the plotted surface declines rapidly (on the scale 
of a few H s ) moving downstream of the impact point. 

2.5. Results for an isothermal equation of state / 
efficient radiative cooling 

Figures 2 through 4 show in greater detail the results 
of the calculation that was illustrated in the left panels 
of Fig. 1. This employed an isothermal equation of state, 
and a ratio of stream to disk scale heights H s /Hd = 2. 
The computational volume extended vertically to 5 disk 
scale heights, and in the x and y directions for 16 disk 
scale heights. We discuss later the differences in simula- 
tions including efficient radiative cooling, but the general 
behavior seen in those runs is the same as in the isother- 
mal case. The computational grid used was 176 x 176 x 55 
cells. We have run simulations including radiative cooling 
at resolutions of up to 192 2 x 60, and adiabatic calcula- 
tions at up to 128 2 x 40. We note, however, that there 
is a sharp vertical density gradient in the disk in the re- 
gion, several scale heights above the midplanc, where the 
stream is overflowing. Higher resolution might thus re- 
veal additional interactions that are not resolved in these 
calculations. 

Fig. 2 shows density contours and velocity vectors in the 
x — y plane, with the stream entering from the bottom in 
each panel, and the disk gas flowing in from the left. In 
the disk midplane (top left panel), the stream is stopped 
rapidly at the edge of the disk by the much denser disk 
material. There is no 'splashing' or ricocheting material. 
Moving upwards, at around 2 disk scale heights above the 
midplanc (top right panel) the densities of the stream and 
disk start to become comparable. Here we resolve the twin 
shocks in the disk and stream gas seen in earlier two dimen- 
sional calculations (Rozcyczka & Schwarzenberg-Czerny 
1987) and predicted analytically (e.g. Livio 1993). We 
also see that the velocity field of the disk gas downstream 
of the impact point has been disturbed by the impact - this 
material now has a significant radial velocity inward. At 
higher elevations the stream rapidly begins to dominate. 
At z = 2.5i/d the stream is already able to penetrate to 
the inner boundary of the computation (at ~ 0.5i? O ut), and 
at z = 3i?d the stream is overflowing almost freely, with 
some modest deflection of the edge facing the oncoming 
disk material. 

Comparing these results with the two dimensional cal- 
culations of Rozcyczka & Schwarzenberg-Czerny (1987), 
which have close to identical resolution in the plane, we 
see the same shock structure and some evidence of unsta- 
ble flow in the mixed stream and disk layer. However our 



3D simulations do not support the idea of a hotspot region 
that is highly extended along the disk rim. In our calcu- 
lations shocks, and hence emission (given the implicit as- 
sumption of rapid cooling), are restricted in the midplanc 
to the part of the disk directly struck by the stream. At 
higher z, the shocks are extended, but primarily along the 
radial direction. The two dimensional nature of earlier cal- 
culations, together with the assumption of a slow fall-off 
of the density at the disk edge, are likely to be responsible 
for these differences. 

Fig. 3 shows density contours in a vertical slice taken 
along the initial stream flow direction. Although the 
stream is converging towards the disk midplane at all radii, 
the interaction between the stream and the disk is concen- 
trated at the disk rim. Inward of the rim, there is a clear 
density minimum between the upper surface of the disk 
and the lower part of the overflowing stream. This is in 
accord with the prediction of Lubow (1989), and it sug- 
gests that the final fate of the stream, not modeled here, 
will be to reimpact the disk near to the point of the bal- 
listic stream's closest approach to the central object. We 
note that the separation between overflowing stream and 
disk seen in this simulation is accentuated by the strong 
flaring of the disk scale height implied by the isothermal 
assumption, cx R 3 ' 2 . However, weak interaction be- 
tween stream and disk interior to the rim is also seen in 
a simulation where we took H& cx i? 9//g , appropriate if the 
disk central temperature scales with the run of effective 
temperature T e ff cx i?~ 3 / 4 (a 'standard' disk). Of course 
stronger interaction would be expected if the disk atmo- 
sphere was much more vertically extended than a gaussian. 

Fig. 4 shows the density on vertical cylindrical surfaces 
at constant radii in the disk (i.e. the co-ordinates are 
R<p,z). Three radii are depicted, at R > i? ut, showing 
just the cross-section of the stream, at the disk rim where 
interaction is strongest, and where the stream is overflow- 
ing the disk interior to the rim. The distortion of the 
stream by the ram pressure of the disk can be seen, and 
it is clear that the main effect of the stream on the disk 
rim downstream of the impact point is to greatly reduce 
the density at high z. The stream acts as a windbreak in 
the flow, and truncates the disk in the vertical direction. 
For these parameters a very small quantity of disk gas is 
deflected over the stream - this phenomenon is stronger 
in simulations where radiative cooling is included and can 
lead to a hot wake of material downstream of the impact 
point along the rim. 

These results are obtained for strictly isothermal flows. 
However, the flow in simulations including explicit radia- 
tive cooling is extremely similar in the regime where the 
cooling is much more rapid than the dynamical timescale. 
As argued earlier, this should be applicable to disks that 
are optically thin at all reasonable accretion rates. 

Fig. 5 compares the results of the isothermal calculation 
with two simulations which include cooling. The cooling 
simulations assume central disk electron densities of 10 14 
and 10 13 particles per cm~ 3 , which correspond to accre- 
tion rates of order 10~ 10 to 10~ n M Q yr _1 . The cooling 
simulations used a 128 2 x 40 grid. A floor temperature 
of 10 4 K was imposed. We plot vertical profiles of density 
and, for the runs with cooling, temperature, at a point 20° 
downstream of the stream impact point just interior to the 
disk rim. 
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The density profiles show that the low density wake be- 
hind the overflowing stream persists in all three runs. In 
the runs incorporating cooling the density at high z can 
be two to three orders of magnitude higher than for the 
isothermal case - this is a consequence of the overspill of 
disk gas above the stream mentioned earlier. The absolute 
density though is very low in all cases. High temperatures 
of 10 5 K and above are seen in a layer above the disk 
surface, this is material that has been shock heated by 
the interaction with the stream and which has passed be- 
neath the overflowing gas. This feature is not adequately 
resolved in these simulations. 

2.6. Results for inefficient cooling 

Fig. 6 summarizes the results for a simulation with an 
adiabatic equation of state and no radiative cooling. The 
grid was again 128 2 x 40 cells, and H s /H d = 2. The re- 
maining initial conditions were as for the isothermal and 
radiative cooling runs, except that the stream injection 
point was moved slightly further from the x = boundary 
in order to reduce effects arising from the finite volume of 
the computational box. 

Qualitatively, the distinction from the isothermal cal- 
culations is that in this simulation the hot shock-heated 
gas expands in all directions and disrupts the coherent 
overflowing stream seen previously. In the disk midplane, 
'splashing' of hot material downstream of the impact point 
is evident, though the velocities of this gas are substan- 
tially less than those of the disk. At z = 3i?d there is 
overflow in a broad fan over the disk surface. Most of 
this material appears to be stream gas that has suffered a 
strong deflection upon reaching the disk edge, though an 
examination of the velocity field also shows disk material 
at high z that is being deflected inwards as a result of the 
stream interaction. 

A slice along the initial stream flow direction now re- 
solves clearly the shocks in the stream and disk gas. The 
overflowing component has a velocity that is predomi- 
nantly in the x — y plane, with a small (a few times c s ) 
vertical component. Further from the disk midplane signif- 
icant vertical velocity away from the disk is seen, suggest- 
ing that with this equation of state a substantially greater 
absorption column will be generated for lincs-of-sight well 
away from the disk plane. As shown later, the actual frac- 
tion of the stream overflowing the disk rim is comparable 
to that in the isothermal case, but the structure of the flow 
is here very distinct. 

2.7. Global effects 

The ZEUS simulations discussed in this paper suggest a 
picture in which efficient radiative cooling leads to a coher- 
ent stream overflowing the disk in a manner very similar to 
that described by Lubow (1989). Conversely, if cooling is 
inefficient, there is a 'splashing' of stream material off the 
disk edge. Similar quantities of material can flow inward 
in this regime, but there is no coherent stream and gas is 
also thrown outwards and upwards. 

In order to examine the qualitative differences in global 
effects between the isothermal and inefficient cooling cases, 
we performed two three-dimensional smooth particle hy- 
drodynamics (SPH) calculations. The results are shown 
in Fig. 7. These calculations have vastly lower resolu- 
tion of the impact region than the ZEUS simulations, but 



cover the whole disk and have been evolved for long enough 
(<~ 10 binary orbits) to allow the disk to relax in the bi- 
nary potential. We show an isothermal calculation (sim- 
ilar to that described in Armitage & Livio 1996), and a 
calculation with a polytropic equation of state (7 = 1.1) 
and system parameters appropriate to the supersoft source 
CAL87 (Callanan & Charles 1989; Gould 1995). 

Similar features are evident in the SPH calculations as 
in the finite difference simulations. The isothermal run 
shows a narrow overflowing stream, while the polytropic 
calculation leads to a much more vertically extended spray 
of material and a bulge of outflowing gas near the impact 
point. Negligible masses of material escape the accreting 
star's Roche lobe in either simulation. Although the choice 
of a value for 7 in such a calculation is essentially arbitrary, 
it is clear from Figure 7 that the trend is such that sys- 
tems where cooling is inefficient should display absorption 
at much higher elevations above the disk plane than those 
well described by an isothermal equation of state. 

3. DEGREE OF STREAM OVERFLOW 

To quantify the degree of stream overflow in the simula- 
tions, we compute the integrated radial mass flux through 
vertical x — z slices which has a radial velocity vr > v cut . 
For w cut < 22c s - the initial inflow velocity of the stream 
gas at the grid boundary - this quantity is conserved in 
the absence of interaction with the disk. Fig. 8 shows the 
radial mass flux as a function of R/R out (where the R for 
each x ~ z slice is evaluated along the line of the initial 
stream flow direction), for a variety of threshold radial ve- 
locities. 

Outside the outer edge of the disk, all the curves coincide 
and are flat, reflecting steady-state undisturbed stream 
flow. For the runs with an isothermal or cooling equa- 
tion of state, the amount of material flowing inward at the 
fastest velocity, > 20c s , drops sharply as the disk rim is 
reached, and the central portions of the gas stream are 
stopped by the denser disk gas. In this case, and as noted 
previously, the main interaction occurs close to the disk 
rim, after which the highest velocity curve remains fairly 
flat at a lower level than initially. This mass flux is com- 
prised of material that has overflowed the disk at large 
enough z to avoid strong interaction with the disk. For 
the illustrated case with Hs/Hd = 2, this is around 5% of 
the stream mass transfer rate. 

From the Figure, it is evident that a significantly larger 
quantity of mass flows inward at velocities that are smaller 
than the ballistic stream velocity, but still highly super- 
sonic. This is either stream material slowed by strong in- 
teraction with the disk, or disk gas entrained by the over- 
flowing stream. For the case shown in Fig. 8, a radial mass 
flux of around 20% of the initial stream mass transfer rate 
is flowing inward with vr > 5c s at the innermost bound- 
ary of the simulation volume, at ~ 0.6 i? ut- We analyze 
the observational signatures of this material more fully in 
the following Section. We also note that close to the rim, 
the total mass of gas flowing inward supersonically exceeds 
the stream mass flux by a factor of ~ 2 — 3, though this 
does not extend far inward. 

For the calculations with an adiabatic equation of state, 
and no cooling, the curves shown in Fig. 8 are different. 
In particular, there is a prominent dip in the curve with 
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the highest velocity cut near the outer rim of the disk, 
corresponding to the position of the shocks seen in Fig. 6. 
However, although the overflowing gas forms much less of a 
coherent stream in this case as compared to the isothermal 
simulations, the total mass moving inward is comparable. 
As with the calculation including cooling, around 10% of 
the stream mass flux overflows the disk at approximately 
ballistic velocities (although the stream is anything but 
ballistic in this limit), with substantially more gas flowing 
inward at lesser velocities. 

For an estimate of the fraction of the stream mass trans- 
fer rate that overflows the disk, we assume that all the 
mass in the stream above some critical z = z CI i t manages 
to overflow the disk, while the stream below that height 
is stopped by the collision with the disk. For a simple 
approach, we take z cr ;t to scale with the height where the 
stream density equals that of the disk (e.g. Livio, Soker & 
Dgani 1986), 

a H d H s ( PdoV^ 2 /o\ 
Zcrit = P —rjo In , (8) 

where we have included a free parameter (5 as a scaling of 
this simple estimate. 

Fig. 9 shows the results of 3 simulations, all incorpo- 
rating radiative cooling, with varying ratios of stream to 
disk scale height. Points are plotted for three radial veloc- 
ity cuts, together with a curve showing the fitting formula 
given in equation (8). As expected, the overflowing mass 
fraction is a strong function of H s /Hd. For these param- 
eters, Hs/H^ = 1.5 implies that only ~ 2 x 10~ 4 of the 
stream passes inward of the disk rim without strong inter- 
action with the disk, although almost two orders of magni- 
tude more material is moving inward at velocities greater 
than 5 times the sound speed. The naive fitting formula, 
with an appropriate choice of (3 (which is however close to 
unity), provides a good approximation to the simulation 
results. 

4. VELOCITY MIXING BETWEEN STREAM AND DISK GAS 

To quantify the observable effects of the stream-disk in- 
teraction on the velocity of the overflowing stream, we 
have computed the expected distribution of radial veloci- 
ties, vr, of gas along lines of sight to the central accreting 
object. We define the elevation angle of the line of sight to 
the disk midplane as 9, and consider viewing the system 
at a phase 18° downstream of the impact point (i.e along a 
diagonal line in Fig. 2 extending from the lower right cor- 
ner towards the central object which is out of the frame 
beyond the upper left corner). We compare the results 
with those expected from a freely flowing stream in the 
absence of the interaction with a disk, obtained by run- 
ning a simulation with the same stream initial conditions 
but no disk. This comparison calculation displays close to 
ballistic flow, but does include stream pressure effects. All 
the calculations in this Section used H s /Hd = 2.5, so that 
a significant fraction of the stream overflows the disk rim. 

Fig. 10a shows the results for an angle of elevation above 
the disk midplane of = 12°. This line of sight in the 
isothermal or cooling simulations probes freely overflow- 
ing gas at high z that has not undergone strong interaction 
with the disk. Consequently, we find that the radial ve- 
locity distribution is almost identical to that which would 



arise from the same ballistic stream if the disk were ab- 
sent. The distribution has a width of ~ 5c s , with a sharp 
cut-off at large negative vr. 

Fig. 10b shows the results closer to the disk midplane, 
at 6 = 9°. This line of sight probes gas which has had 
significant interaction with the disk. The interaction leads 
to reduced (less negative) radial velocities in the column 
of gas, corresponding to the ram pressure induced bending 
of the inner edge of the stream seen in Fig. 1, and a nar- 
rower distribution of vr. The deviations from the ballistic 
prediction for the line of sight radial velocity are here at 
the - 10% level. 

In general we find that for these relatively large angles 
above the disk midplane, the column of gas at large neg- 
ative vr greatly exceeds that of the undisturbed disk at 
vr w 0. These velocity differences with the ballistic ex- 
pectation are therefore potentially observable as absorp- 
tion against the central object. For lines of sight that are 
closer to the disk midplane, we find gas with a range of 
radial velocities that become steadily closer to vr = as 
9 decreases. However, there is always a substantial com- 
ponent at vr — which arises from the undisturbed disk 
interior to the overflowing stream. This much lower veloc- 
ity gas would therefore only be observable in emission. 

Finally, we show in Fig. 10c the result for an adiabatic 
equation of state. As expected from the flow pattern seen 
in Fig. 6, the adiabatic equation of state leads both to 
a greater column of absorbing gas at high angles above 
the disk midplane, and a much larger spread of velocities 
of that gas. As compared to the isothermal simulations, 
the mean velocity of the gas along the line of sight is ap- 
proximately halved, and there is a very broad spread of 
velocities extending even to positive vr. 

5. COMPARISON WITH OBSERVATIONS 

Three classes of observations potentially probe the in- 
teractions modeled here: X-ray absorption 'dips' in nearly 
edge-on systems, which may track the column of over- 
flowing gas; observations of the hotspot region where the 
shocks occur; and doppler tomography analyses which may 
be able to distinguish the velocities predicted to arise as 
a consequence of the stream impact. We briefly discuss 
these in turn. 

Dips are seen in the UV and X-ray light curves of nearly 
edge-on low-mass X-ray binaries (LMXBs - e.g White, Na- 
gase & Parmar 1995), intermediate polars, and dwarf no- 
vae (Long et al. 1996; Szkody et al. 1996). In the case of 
LMXBs, these have been interpreted as arising from the 
stream-disk interaction, either in the form of an overflow- 
ing component (Frank, King & Lasota 1987), or as a bulge 
excited on the rim of the disk (Hellier & Mason 1989; Ma- 
son 1989). Similar features seen in CVs also seem to be 
consistent with this interpretation. In particular, in Z Cha 
the dips in a UV lightcurvc were seen to become shallower 
during the course of a superout burst, which would be con- 
sistent with a declining mass transfer rate through the 
event (Harlaftis et al. 1992). For the dips observed in the 
intermediate polars FO Aqr, EX Hya, BG CMi and AO 
Psc, and in the dwarf nova U Gem in X-ray light curves 
(Hellier, Garlick & Mason 1993), the dips are observed to 
become shallower with decreasing orbital period of the sys- 
tems. This is consistent with a statistically smaller mass 
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transfer rate in shorter period systems (e.g. Warner 1987; 
Patterson 1984), and an interpretation of the dips as aris- 
ing from the stream-disk impact. It is also possibly fortu- 
itous, as both the simulations and common sense suggest 
that the absorbing column should be a strong function of 
the system inclination. Moreover, these observations do 
not discriminate between the overflow and bulge models 
for the absorbing column, and there remains the obser- 
vation of dips at a plethora of phases other that ~ 0.8 
(e.g. Mason 1989; Thorstensen et al. 1991), where the im- 
pact occurs. Some of these may arise at the point where 
the stream converges and reimpacts the disk, which oc- 
curs at a binary phase of ~ 0.6 for a variety of mass ratios 
(Lubow 1989), or be generated due to an eccentric disk in 
low mass ratio systems. The latter possibility is suggested 
by the observations of Harlaftis et al., and supported by 
theoretical calculations (Whitehurst 1988a, 1988b; Lubow 
1991, Armitage & Livio 1996; Murray 1996). 

We also note that models for the structure of the 
accretion disk rim in supcrsoft X-ray sources require 
significant 'splashing' at the stream-disk impact point 
(Meyer-Hofmeister, Schandl & Meyer 1997). Since the 
accretion rates in these systems are extremely high (~ 
10~ 7 Mgyr" 1 , e.g. van den Hcuvel et al. 1992), this 
requirement is consistent with the results presented here. 
The effects of strong irradiation from the accreting object 
are an additional important complication in both these 
systems and in low-mass X-ray binaries (see, e.g the sim- 
ulations reported by Blondin 1997). 

The simulations reported in this paper suggest that 
emission from the pair of shocks generated by the impact 
of the stream with the disk occurs over a smaller azimuthal 
extent than was implied by earlier calculations restricted 
to two dimensions. This is in general agreement with ob- 
servations, for example of Z Cha, where Wood et al. (1986) 
obtain a hotspot extent of ~ 4° in azimuth, and IP Peg, 
where Wood & Crawford (1986) find an extent of ~ 3 — 7°. 
The radially extended shock structure seen in the simula- 
tions is consistent with an apparently greater azimuthal 
extent of the hotspot when seen in ingress as compared to 
at egress, this is seen in IP Peg (Wood & Crawford 1986). 

Optical spectroscopy provides a further probe of the ab- 
sorption and emission expected to arise from an overflow- 
ing stream. For example, Hellier (1996) reports observa- 
tions of the nova-like variable V1315 Aql. In this system, 
absorption that has been attributed to stream overflow is 
seen at phases close to the point where the stream impacts 
the disk. Further along the stream trajectory, emission is 
seen at velocities lower than the free-fall prediction. This 
would appear to be consistent with an approximately bal- 
listic overflow model. However doppler mapping of an- 
other nova-like variable, SW Sex, shows emission that, if 
attributed to the stream-disk impact, seems more consis- 
tent with a bulge along the disk rim (Dhillon, Marsh & 
Jones 1997). This is also the interpretation advanced for 
doppler maps of the dwarf nova WZ Sge (Spruit & Ruttcn 
1997), where the accretion rate is much lower. 

Clearly, if the interpretation of all of these observations 
is correct, they do not fit nicely into the simple picture 
presented earlier, where the accretion rate determines the 
geometry of the stream overflow. However it is at least 
encouraging that observations, for example of OY Car 
(Billington et al. 1996; see also Bruch, Beele & Baptista 



1996), may be capable of constraining the radial location 
of absorption features. This, together with theoretical pre- 
dictions of the velocity fields arising from the stream-disk 
interaction, should help to disentangle the structure of the 
emission and absorption regions generated by the impact. 

6. SUMMARY AND CONCLUSIONS 

In this paper, we have discussed high resolution simu- 
lations of the stream-disk interaction in binary systems. 
We find as our main result that the qualitative behav- 
ior of the flow depends on the efficacy of cooling in the 
shock heated gas created by the impact. If cooling is effi- 
cient, specifically if the cooling time is much shorter than 
the dynamical timescale in the disk, then the up- 
per reaches of the stream freely overflow the disk. This 
approximately ballistic flow is as described by earlier cal- 
culations (e.g. Frank, King & Lasota 1987; Lubow 1989), 
and is expected to terminate in a secondary hotspot where 
the stream reimpacts the disk close to its closest approach 
to the accreting object. We see no evidence for hydro- 
dynamic instabilities leading to rapid entrainment of the 
entire stream as it flows over the disk. Additionally we find 
significantly larger quantities of mass that undergo strong 
interaction with the disk (or, are entrained disk gas), and 
flow inward with intermediate velocities. This gas may 
dominate the observed apparent velocity of the overflow- 
ing material at intermediate angles of disk inclination, and 
is likely to be visible in emission. We have argued that this 
efficient cooling mode of stream overflow is appropriate to 
CVs at low accretion rates, where the hotspot emission is 
emanating from regions that are not optically thick. 

Conversely, if cooling is inefficient, then the interaction 
leads to a flow better envisaged as an explosion at the 
point where the stream impacts the disk (e.g. Livio, Soker 
& Dgani 1986). Material expands away from this point in 
all directions, and there is no coherent stream flow over 
the disk. Similar fractions of the incident gas stream flow 
inward in either regime, but in the case of inefficient cool- 
ing the impact leads to a more vertically extended 'halo' 
above and below the disk, and a bulge extending along the 
disk rim. The velocity of gas viewed along a line-of-sight 
to the accreting star displays a much greater dispersion in 
the case of inefficient cooling, and has a mean velocity that 
has a smaller magnitude than if cooling is efficient and the 
stream overflows ballistically. 

Aside from the strength of cooling, the most important 
parameter in determining the degree of stream overflow 
is likely to be the ratio of stream to disk scale heights at 
the disk edge. We have quantified this for the case of ef- 
ficient cooling. Assuming, as we have here, that both the 
stream and the disk have vertical density profiles that can 
be described as gaussians, then H s /H^ ^ 2 is required for 
significant stream overflow at ballistic velocities to occur 
(at lower ratios significant supersonic inflow at lesser ve- 
locities may still be produced). We find that the numerical 
results can be reasonably fit by a simple expression that 
assumes that the overflowing mass fraction scales with the 
fraction of the stream above the height where the disk and 
stream densities are equal. 

The simulations reported here have used very simple 
isothermal vertical disk profiles, and this is probably the 
largest remaining uncertainty in the calculations. A den- 
sity profile that declines markedly less steeply than the 
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isothermal one would lead to stronger interaction between 
the stream and the disk at several scale heights above the 
midplanc. The observations discussed above suggest that 
overflow can occur in CVs, but in extreme cases such a 
density profile could inhibit overflow entirely Even with 
weaker interactions, the influence of the stream on a disk 
wind or on large-scale magnetic structures anchored to the 
disk surface might extend the non-axisymmetric effects of 
the stream to larger elevations above the disk surface than 



the purely hydrodynamic effects studied in this paper. 
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Fig. 1. — Isodcnsity surfaces from the isothermal (left column) and adiabatic (right column) simulations. The upper panels show the top 
view, with the stream flowing left to right and the disk flowing downwards; the lower panels show side views in which the disk material is 
moving out of the plane of the paper. The isothermal calculation is rendered at a density surface that is 10 — 3 of the central disk density. 
The adiabatic run is rendered at a density surface that is 10~ 2 ' 75 of the central disk density. This figure available as a jpeg in the astro-ph 
version. Interactive viewing of the simulations at http : //www. cita.utoronto . ca/~ armitage/hydro_abs .html. 




Fig. 2. — Density contours in the x — y plane for a calculation with an isothermal equation of state. The computational grid was 176 X 176 X 55 
cells. H B /H,i = 2. The panels depict slices in the disk midplanc (z = 0), and at 2, 2.5 and 3 disk scale heights above the plane. The stream 
enters from the bottom in each panel, the disk flow is left to right. Contour levels are at Alogp = 0.3, between logp = and logp = —12. 
Arrows depict the direction and relative magnitude of the velocity field in each plane. The radius as a fraction of the disk radius Rout is also 
shown. 
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Fig. 3. — Density contours in a vertical slice along the initial stream flow direction for the isothermal calculation. Contour levels are at 
Alogp = 0.3. 
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Fig. 4. — Density contours on cylindrical surfaces at constant radius from the accreting object (R(j>,z surfaces). From top downwards, the 
surfaces are at radii of 1.2R ou t, 1.0-Rout and 0.75-R O ut- Contour levels arc at A logp = 0.3, between logp = and logp = —9. The rr-axis is 
labelled with the angular distance relative to where the stream meets the disk edge — in this figure positive A<p implies downstream of the 
impact point. 
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Fig. 5. — Illustration of the influence of the cooling rate on the vertical density and temperature profile downstream of the impact point 
on the disk rim (at R = 0.9R ou t, and 20° downstream of the stream impact position). The solid line shows the result for an isothermal 
equation of state, the short dashed and long dashed lines for calculations including radiative cooling with central disk densities of f0 14 cm~ 3 
and 10 13 cm -3 respectively. The dotted line in the density profile depicts the hydrostatic fall-off of density with height. 
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Fig. 6. — Summary of results for an adiabatic equation of state, with initial conditions corresponding to the isothermal calculation. The 
computational grid was 128 X 128 X 40 cells. The upper panels show density contours in the x — y plane at z = and at 3 disk scale heights, 
the lower panel a vertical slice along the initial stream flow direction. Contour levels are at Alogp = 0.2. Arrows depict the direction and 
relative magnitude of the velocity field in each plane. Note that the vertical slice is not in the radial direction, this is why there are non-zero 
velocities even within the disk material. 



Fig. 7. — Results from SPH calculations of the entire disk with varying polytropic equations of state. Plotted is a side view of the system 
with the stream entering from the left. Upper panel is for an isothermal equation of state, lower panel from a calculation with 7 = 1.1. 
Omitted for reasons of space. Armitage & Livio (1996) shows a similar figure. 
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Fig. 9. — Overflowing mass as a fraction / of the stream mass flux, as a function of the ratio of scale heights of the stream to the disk 
-Hstrcam /-Hdisk- The symbols are from the numerical simulations; crosses for mass inflow at vr > 5c s , squares for vr > 10c s , and triangles 
for vr > 20c s . The curves represent the scaling function described in the text, with (3 = 0.85, 1.1 and 1.5 respectively. 
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Fig. 10. — The radial velocity distribution of gas along lines of sight to the central accreting object. All panels are for a line of sight situated 
18° downstream of the impact point, (a) Elevation angle above the midplanc 9 = 12°. Solid line is the result from a simulation including 
cooling, dashed line is the contribution from a ballistic stream with no disk interaction, (b) As in the previous panel except with 9 = 9°. (c) 
Comparison of an adiabatic (long dashed line) simulation with a simulation including efficient cooling (solid line), for 9 = 12° (note different 
scale on the x axis). 
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